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Parity violation in the Zff coupling can be measured via the observables Af = 
2vfdf/(v'f + ai), where Vf and af represent the vector and axial vector couplings to fermion 
/. In particular, for f = b, Af, is largely independent of propagator effects that modify the 
effective weak mixing angle, and thus provides an unambiguous test of the Standard Model. 

The Born-level differential cross section for the process e + e _ — ► Z° — > // is 

da f /dz<x(l- A e P e )(l + z 2 ) + 2A f (A e - P e )z , (1) 

where P e is the e~ beam longitudinal polarization (P e > for right-handed (R) polarization) 
and z is the cosine of the polar angle of the outgoing fermion with respect to the incident 
electron. The ability to modulate the sign of P e allows the final-state quark coupling Af 
to be extracted independently of A e from a fit to the differential cross section. Thus, the 
measurements of Af described here are unique, and complementary to other electroweak 
measurements performed at the Z° pole 0. 

This Letter reports the results of the 1996-98 SLD lepton tag analysis, for which identified 
electrons and muons were used to tag the flavor of the underlying heavy quark. The data 
sample used in this analysis is roughly three times larger than that of previously reported 
results 0. Further statistical and systematic advantage is provided by improvements to the 
data analysis which take advantage of the precise information provided by the new vertex 
detector (VXD3) that was installed just prior to the 1996 data run. 

The Stanford Linear Collider (SLC) and its operation with a polarized electron beam 
have been described elsewhere Q. During the 1996-98 run, the SLC Large Detector (SLD) || 
recorded an integrated luminosity of 14.0 pb _1 at a mean center of mass energy of 91.24 
GeV, with a luminosity-weighted electron beam polarization of \P e \ = 0.7336 ± 0.0038 ||. 

Charged particle tracks are reconstructed in the Central Drift Chamber (CDC) and the 
CCD-based vertex detector in a uniform axial magnetic field of 0.6T. For the 1996-98 data, 
the combined CDC and VXD3 impact parameter resolution in the transverse (longitudinal) 
direction with respect to the beam is 7.7 (9.6) /im at high momentum, and 34 (34) /im at 
p±V sin 9 = 1 GeV/c 2 , where p± and 6 are the momentum transverse to and angle relative to 
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the electron beam direction. The Liquid Argon Calorimeter (LAC) measures the energy and 
shower profile of charged and neutral particles with an electromagnetic energy resolution 
of cfe/E = 15%/ \J E{GeV) and is used in the electron identification. The Warm Iron 
Calorimeter (WIC) detects charged particles that penetrate the 3.5 interaction lengths of the 
LAC and magnet coil. The Cherenkov Ring Imaging Detector (CRID) measures the velocity 
of charged tracks in the region | cos#| < 0.68 using the number and angle of Cherenkov 
photons emitted in liquid and gaseous radiators; electrons are well separated from pions in 
the region between 2 and 5 GeV/c, while pion (kaon) rejection reduces backgrounds to the 
muon sample in the region 2 < p < 5 (2 < p < 15) GeV/c. 

The axis of the jet nearest in angle to the lepton candidate is used to approximate z, the 
cosine of the polar angle of the underlying quark. Jets are formed from calorimeter energy 
clusters (including any associated with the lepton candidate) using the JADE algorithm 
with parameter y cut = 0.005. The analyses presented here make substantial use of 
'secondary' decay vertices which are displaced from the primary interaction point, identified 
via the ZVTOP topological vertexing algorithm ||, as well as the invariant mass of the 
tracks comprising the secondary vertex ('vertex mass'), corrected to account for unmeasured 
neutral particles 0. 

The selection of electron and muon candidates with p > 2 GeV/c in hadronic Z° de- 
cays has been described previously @. Electrons are identified with both LAC and CRID 
information for CDC tracks in the angular range | cos#| < 0.72. Electrons from photon con- 
versions are recognized and removed with 73% efficiency WIC information is also included 
for muons, providing an essential measurement of their penetration. Muons are identified 
in the angular region | cos#| < 0.70, although the identification efficiency falls rapidly for 
| cos#| > 0.60 due to the limited angular coverage of the WIC. To reduce backgrounds from 
misidentification, the 29% of events containing electron candidates that had no reconstructed 
secondary vertices were removed from the sample, precluding the use of the electron sample 
for the measurement of A c . 

For p > 2 GeV/c, Monte Carlo (MC) studies indicate efficiencies (purities) of 64% 
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(64%) and 81% (68%) for the electron and muon samples, respectively, where the remaining 
electrons from photon conversion account for 5% of the 12862 electron candidates. In the 
case of the muon sample (21199 candidates), the background is due both to misidentification 
(8% of muon candidates) and to real muons from light hadron decays (25%). In both cases, 
the MC simulation has been verified with a control sample of pions from Kg — > tt + tt~ decays. 
The fraction of such pions misidentified as electrons is (1.02 ± 0.06)%, consistent with the 
MC expectation of (1.06 ± 0.03)%. For muons, the measured pion misidentification fraction 
is (0.342 ± 0.028)%, somewhat higher than the MC expectation of (0.279 ± 0.012)%. This 
difference has been accounted for by raising the background level in the maximum likelihood 
fit to the muon sample by (20 ± 10)% of itself. 

The sample of events containing identified leptons is composed of the following 
event types (charge conjugates implied): Z° — > bb,b — > I ('&/'); Z° — > bb, b — > c — > I ('fee/'); 
Z° — > bb,b — > c — > I ( l bcV); Z° — > cc, c — > I ('c/'); and background from light hadron and vec- 
tor meson decays, photon conversions, and misidentified hadrons ('6/c'). 

Identification of electron candidate event types is based on the values of eight discrim- 
inating variables JTTJ: track momentum (p), momentum transverse to the nearest jet (p t ), 



the estimate of the underlying B hadron boost |TT| and when available, same hemisphere 
secondary vertex mass, opposite hemisphere vertex mass, same hemisphere vertex momen- 
tum resultant, same hemisphere vertex significance (separation D between the interaction 
point and secondary vertex, divided by its uncertainty), and L/D (where L is the distance 
from the interaction point to the point on the secondary vertex trajectory closest to the 
electron candidate trajectory). These variables are used as inputs to an Artificial Neural 
Network with three output nodes Nu, N^i, and N c i, optimized for the bl, bcl + bcl, and cl 
signals, respectively. Event type probabilities are estimated according to the composition of 
MC electron candidate events with similar output node values. The measured and simulated 
distributions of the three output node variables are compared in Figure |1|. 

The Neural Network is trained on the SLD MC sample of hadronic Z° decays, generated 
with JETSET 7.4 [12|. Semileptonic decays of B mesons are generated according to the 



ISGW formalism Jl3| with a 23% D** fraction, while semileptonic decays of D mesons are 
simulated according to branching ratios reported by the Particle Data Group ||14|| . Exper- 
imental constraints are provided by the B — > I and B — > D inclusive momentum spectra 



measured by the CLEO collaboration [|I5||L6|| and the D — > I momentum spectrum measured 



by the DELCO collaboration 17]. The detailed simulation of the SLD detector response 
has been realized using GEANT ||18|| . 

Muon candidate event type probabilities are estimated according to the composition of 



MC muon candidate events with similar values of the following discriminating variables |19 
p, pt, and, when available, L/D and M max , the largest of the secondary vertex invariant 
masses. The measured and simulated distributions of these variables are compared in Fig- 
ure ||. 

A maximum likelihood analysis of all selected hadronic Z° events containing lepton 
candidates is used to determine A b and A c . The likelihood function contains the following 
probability term for each lepton, with measured charge sign Q: 

P (P e , z; A b ) ex {(1 + z 2 ){l - A e P e ) - 2Q(A e - P e ) [(/ w (l - 2x b ) - f M (l - 2 Xb5 ) 

+ f m (l - 2x 65 ))(l - A b QCD (z))A b - f- d (l - A c QCD (z))A c + f bk A bk ] z] . (2) 

The lepton source fractions fu, f b5l , f b5 i, f 5 i, and f bk are functions of the three neural 
net output node values (electron candidates) or the four discriminating variables (muon 
candidates). For the fit to muon candidates, both A b and A c are left as free parameters, 
whereas A c is fixed to its SM value (see Table I) for the fit to electron candidates. 

Correction factors (l—2x x ), where \ x is the mixed fraction for lepton source x, are applied 
to b-quark lepton sources to account for asymmetry dilution due to B°B° mixing. The value 
of Xb is taken from LEP measurements of the average mixing in semileptonic B decays [|TJ, 
but must be corrected to take into account selection and fitting bias, including that due 
to the enhanced likelihood for bcl cascade leptons to have come from a B meson which has 
mixed |2"0fl . For the electron sample, MC studies indicated that the mixing probability Xb for 
bl decays was independent of of the value of the NN output parameters, but was increased 



by a relative 1.7% overall by the bias of the vertex requirement towards the selection of B° 
over B ± decays. For the muon sample, the effective values of xis an d Xbc were evaluated on 
an event-by-event basis, based on MC events with values of the muon-sample discriminating 
variables close to those of the given data event. 

The asymmetry in the background is parameterized as a function of p and p t . For the 
electron sample, the parameterization is determined from tracks in the data not identified 
as leptons. For the muon sample, MC studies indicated a substantial difference between 
the true background asymmetry and that of non-leptonic tracks, and so the background 
asymmetry parameterization was determined directly from the MC simulation. 

A ^-dependent correction factor (1 — Aq CD (z)) is included in the likelihood function to 
incorporate the effects of gluon radiation. Calculation of the quantity Aq CD (z) has been 



performed by several groups pl| . For an unbiased sample of bb or cc events with \z\ < 0.7, 
correcting for this effect increases the measured asymmetry by ~ 3% overall. However, a 
MC simulation of the analysis chain indicates that biases which favor qq events over qqg 
events mitigate the effects of leading order gluon radiation by about 30%. Effects due to 
gluon splitting to bb and cc have been estimated by rescaling the JETSET simulation to 



world average gluon splitting measurements [Q . Additional radiative effects, such as those 
due to initial-state radiation and 7/Z interference, lead to a further correction of —0.2% 
(-0.1%) on the value of A b (A c ). 

A list of systematic errors is shown in Table |. The purity of the separation of Z° — > bb 
and Z° — > cc events via secondary vertex information introduces an uncertainty dominated 
by the efficiency of charged track reconstruction, which has been constrained by reweighting 
MC tracks by the ratio of the number of tracks in data and MC as a function of p and pt- 
The ability of the L /D variable to discriminate between bl and bcl decays is sensitive to the 
fraction of B — > DD decays, which has been constrained from SLD data [[19]. 

For the 1996-98 muon sample, we find that A b = 0.938 ± 0.044(stat.) ± 0.024(syst.) and 
A c = 0.560 ± 0.063(stat.) ± 0.064(syst.), with a statistical correlation coefficient of 0.108. 
For the corresponding electron sample, we find A b = 0.896 ± 0.050(stat.) ± 0.028(syst.). 
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Combined with the result of , we find overall SLD average results via semileptonic B and 
D hadron decay of 

A b = 0.919 ± 0.030(stat.) ± 0.024(syst.) 
A c = 0.583 ± 0.055(stat.) ± 0.055(syst.). 

In conclusion, we have directly measured the extent of parity violation in the coupling of 
Z° bosons to b and c quarks using identified charged leptons from semileptonic decays. The 
results presented here take advantage of an additional sample of 400000 Z° decays, and em- 
ploy a new method of signal source separation, resulting in substantial increases in precision 
relative to previous measurements f2|. These results are in agreement with the Standard 
Model predictions Af, = 0.935 and A c = 0.667, which are insensitive to uncertainties in 
Standard Model parameters such as the strong and electromagnetic coupling strengths, and 
the top quark and Higgs boson masses. 

We thank the staff of the SLAC accelerator department for their outstanding efforts 
on our behalf. This work was supported by the U.S. Department of Energy and National 
Science Foundation, the UK Particle Physics and Astronomy Research Council, the Istituto 
Nazionale di Fisica Nucleare of Italy and the Japan-US Cooperative Research Project on 
High Energy Physics. 
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